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A new class of exact solutions of hydrodynamic equations for an incompressible fluid (gas) 
at the presence of a bulk sink and uprising vertical flows of matter is considered. The 
acceleration of the rotation velocity of classical non-stationary vortices is conditioned by 
the joint action of the convective and Coriolis hydrodynamic forces (accelerations), which 
appear due to the converging radial flows of the matter in the region of a bulk sink. It 
is shown that there exist velocity proflles that nullify viscous terms in the Navier-Stokes 
equations and represent a vortex structure with a " rigid-body" rotation of its core and 
converging radial flows. 

The concept of non-stationary vortices in open systems is applied to description of 
origination of power atmospheric vortices (whirlwinds, tornados, and typhoons). In the 
classical hydrodynamics a favorable condition for the origination and existence of such 
vortices is the exact nulliflcation of the terms, which describe kinematic viscosity of 
an incompressible fluid. Such flows have the minimal rate of energy dissipation that 
corresponds to the " minimum entropy production principle" , and therefore may relatively 
easily appear in favorable natural conditions. 



1. Introduction 

Hydrodynamical vortices and vortex-like flows are extremely widespread in nature. 
This is evidently proven by constant occurrence of powerful atmospheric vortices — cy- 
clones, hurricanes, typhoons, tornados (see Nalivkin (1969)[ Bengsson & Lighthill (1982) ) 



sandstorms, oceanic vortical flows — the so-called "rings", twirled streams (whirlpools) 
on rivers, turbulent vortices in wake swirls of ships, etc. An example of a huge station- 
ary vortical structure is the big "Red Spot" on the surface of Jupiter, which has been 
observed by astronomers for several hundred of years. 

One may ask a natural question: why such hydrodynamic vortex structures can exist for 
a long time in liquids and gases despite the finite viscosity of the medium? The answer to 



this question is partly given by the Rankine model (see Kundu (1990) ) for a vortex in an 
incompressible viscous fluid where its azimuthal velocity profile is set to be linear in the 
radius r of a " rigid-body" rotation v^p (r) = cur with the angular velocity tu inside a certain 
cylindrical region r < i?o of a radius Rq, and have the "differential" rotation Vip{r) = 
ivRl/r in the external region r > Rq. In this case, the terms with a bulk viscosity in 
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the Navier-Stokes equations for an incompressible liquid (see Landau & Lifshits (1987) I 
are precisely nullified in the cylindrical system of coordinates. At the first sight, this 
corresponds to the non-dissipative vortical rotation of a liquid with a constant value of 
Lv = const and with a sharp break of a profile v^{r) at the point r ~ Rq. 



However, the radial dependence of the velocity v^(r) of real hydrodynamical vortices 
cannot have any singular points (in particular, jumps or breaks that correspond to jumps 
of the first derivative) and also must be a continuous analytical function of the radius r. 
Therefore, in a neighborhood of the point r = Rq, the function v^{r) should be smooth, 
and there is a finite viscosity in this area, which results in dissipation of kinetic energy of 
a vortex, its spreading in space and decaying in time. Stationary or accelerating vortices 
can exist only due to inflow of energy from the external environment. At this, from all 
possible hydrodynamical flows and vortical movements in an incompressible liquid, easily 
arise, longer exist and dominate over others only those ones, whose viscosity effects are 
almost absent and energy dissipation plays an insignificant role. Such hydrodynamical 
flows and vortices correspond to the " minimum of entropy production principle" and are 
frequently enough observed in the nature. 



In the present paper we will consider several examples of origination of non-stationary 
vortical structures in open multicomponent systems (solutions, homogeneous and het- 
erogeneous mixtures), in which due to chemical or phase transformations a bulk sink 
(convergence) of the matter exists when one or several components drop out from the 
general collective hydrodynamical motion of the system. In this case, it is supposed that 
the chemical compound of the substance, as well as its density p = const, are constant 
in time and almost homogeneous in space due to dynamical and chemical balance with 
the surrounding environment. 



We will show, that the growth of the rotation velocity of classical non-stationary vor- 
tices is caused by the joint action of convective and Coriolis hydrodynamic forces (accel- 
erations) that arise due to converging radial flows of the substance flowing to the area of 
a bulk sink. The law of the velocity growth in time can be exponential, or can correspond 
to a nonlinear "explosive" regime of instability, when velocities formally reach infinite 
values in a finite time interval (under the condition of unlimited inflow of the substance 
from the environment). 



A limit of the growth of velocity and kinetic energy of a vortex in a normal liquid (gas) 
in real conditions is caused by several factors: the friction between the liquid (or gas) and 
fixed solid surfaces; the development of turbulent dissipative processes, in particular, on 
the border of the vortex core at r = Rq, where a tangential jump of the velocity arises, 
and a connected with it small-scale instability of surface waves; the decrease of inflow of a 
substance from the outside, etc. We shall also note that at approaching of hydrodynamical 
velocity to the speed of sound, the compressibility of a liquid must be taken into account, 
which in turn accounts for dissipation due to the volumetric viscosity. All these factors will 
be considered and estimated below, but before that a simple mechanism of origination of a 
vortex in an incompressible one-component liquid (gas) due to existence of parallel to the 
vortex axis hydrodynamic flows, whose velocity depends on the longitudinal coordinate 
z, will be considered. Such a mechanism may explain the formation of a ftinnel in water 
or a sand-spout in a desert due to accelerated falling of the water in the hole or rising of 
the hot air in the gravitation field of the Earth. 
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2. One-component liquid with accelerated flows: the mechanism of 
formation of a funnel and a windspout 

One of the most interesting paradoxes in hydrodynamics is the so-called "funnel ef- 



fect" (see Kundu (1990) Sedov (1997)). It is usually assumed, that this effect is caused 
by conservation laws of the angular momentum of an incompressible fluid (gas) inside 
a given contour, and is accompanied by the accelerated rotation of a vortex at con- 
centration of the vorticity of a flow ui — rotv due to the narrowing of the channel. 
Another approach to the problem of a funnel formation lies in the assumption (see 



Goldshtik, Shtern & Yavorsky (1989)) about origination of the angular momentum at 
zero initial vorticity as a result of instability of cylindrically-symmetric flow in a liquid 
(the flooded jet) in relation to axially-asymmetric left- and right-spiral perturbations 
with ejection of rotation of a certain sign to the inflnity at the expense of a flow (con- 
vective instability) and accumulation of a rotation motion of the other sign (absolute 
instability). 

Let us show, that there exists one additional simple mechanism of a vortex formation 
in an incompressible liquid (gas), which is present in a gravitational field and includes 
vertical ascending or descending flows, whose velocities depend on the coordinate z (along 
the vertical axis of the vortex) . 

Consider the Navier-Stokes equations for an axially-symmetric motion of an incom- 



pressible viscous fluid (gas) in cylindrical coordinates (see Landau & Lifshits (1987)): 



dVr dVr 1 dP / d'^Vr 1 dVr 



dt dr r p dr \ dr"^ r dr 

'd^v^ ^ 1 dv^ 
r dr 
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p dz 



(2.1) 

(2.2) 
(2.3) 



dvz dvz dvz 1 dP I d Vz 1 dvz d Vz 

+'"r-^ +Vz^ = -T. g + V\ -TTY ^ + ^"l" 

dt dr dz pdz \ dr^ r dr dz-^ 

where v^-^v^ and Vz are radial, azimuthal and axial components of the hydrodynamical 
velocity v, P and p are the pressure and the density of a liquid (gas), v — r\l p is the 
coefficient of kinematic viscosity, and g is the acceleration due to gravity, which is directed 
to the opposite direction of the axis z. 

Equations (|2.1[) - (|2.3[) are completed with the continuity equation 



dVr Vr dVz 

divv= -^ + - + -^ = 0. 2.4 

dr r dz 

Let us notice, that in equations (|2.ip and (12. 2p . dependencies of Vr and on z are not 
taken into account only for simplification. 

In the case of a plane vortical rotation, when Vr = Vz = 0, equations (|2.ip and (|2.2p 
become: 



1 dP 



r 



p dr 



dv^ _ ^ / dhj^ ^ Idv^ 
dt V 9?'^ r dr r'^ 



(2.5) 



(2.6) 



If the radial dependence of the azimuthal velocity Wi^(f), according to the Rankine 



model of a vortex (see Kundu (1990)), is chosen to be 



Mr)-< ""^p"' (2.7) 
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then the right hand side of equation (12. 6|) is identicaUy equal to zero, which corresponds to 
the rotation of an incompressible liquid (gas) with a constant angular velocity uj = const. 

At this, the distribution of hydrodynamical pressure P{r), according to is in the 

so-called 'cyclostrophic' rotation regime and has the form: 

where Pq — Poo — puPRq is the pressure on the vortex axis, and Poo is the pressure at 
large distances (i.e. at r — ^ oo). 

In the Rankine model, the radius of the vortex core i?o is not determined, but if in a 
liquid (gas) there exists a cylindrically-symmetric flow (the flooded jet) with a velocity 
Vz, then its radius will determine the size i?o in equation (|2.7p . Let us assume, that the 
flow velocity depends linearly on the coordinate z and does not depend on r in the region 
r ^ i?o, i.e. we have the profile: 

In this case the continuity equation (|2.4p is satisfied for the following radial dependence 
of velocity Vr, which is continuous aX r — : 



Vr{r) ^{ / R2 (2.10) 




Let us notice, that the above mentioned structure of the velocities (|2.7p , p.9|) and (I2.10p 
identically nullify the viscous terms in equations (|2.1|) - ()2.3|) . At the same time, the di- 
agonal components of the viscous stress tensor are different from zero for these profiles, 
which results in the following expression for the rate of changing of the kinetic energy of 
a vortex due to dissipation (per unit of length of a vortex along the axis z): 



K (2.11) 



Substituting expressions (12. 7|) and ()2.10p in equation ()2.2p . we get 

dio ( auj, r ^ i?o, 



dt ~ \ 0, r> Ro. (^-^^^ 

Equation ()2.12|) is a result of that fact that in the inner region (r ^ Rq) the convertive 
and Coriolis forces at 7^ are added up, while in the external region (r > Rq) they 
mutually compensate each other. 

dck 

If the parameter a is constant in time (a const, — = 0) and positive a > 0, then 

dt 

from equation (|2.12p it follows that inside the region r ^ i?o the angular velocity of a 
vortex grows in time according to the exponential law 

w(t) = w(0)e"*, (2.13) 

provided that a nonzero initial vorticity uj{0) ^ exists in the liquid (gas), whereas in 
the external region lu = uj{0) = const . Thus, on the border of a vortex core r = Rq, a, 
jump of the azimuthal velocity that exponentially grows in time arises (see section |4?2|) . 

Let us emphasize, that the energy dissipation p. lip of the non-stationary vortical 
motion does not depend on time t, and is determined only by the initial vorticity 2a;(0). 
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Thus, the dissipation remains small, despite the fast increase of the angular velocity of 
a "rigid-body" rotation of the vortex core. This means that non-stationary vortices are 
not suppressed by the dissipation at least at an initial stage of their developments. 

On the other hand, substituting expressions (|2.7p . (|2.9p and (|2.10p in equations (|2.ip 
and (j2.3p . we get the following equations for the determination of the hydrodynamical 
pressure: 

a' it) 



dP 
dr 

dP 
dz 



pr 

pi4 



u\t) 



HO) 



1 da 

' 2~dt 
''(t) Idar^ 

~ ^ 2~dtR^ 



g + avzo + z 



da 
Itt 



-P9, 



r > Ro, 

r ^ Rq, 
r > Rq. 



(2.14) 



(2.15) 



Let us notice, that equations (|2.14p and (|2.15p indicate on the existence and increasing 
in time of a jump of the first derivatives of the pressure dP/dr and dP/dz on the surface 
of the vortex core r ^ Rq. From equation (|2.15p . a possibility of existence of a non- 
stationary solution with ^ 7^ follows, namely: 



dP 

^2(0 = 0, —±pg = Q. 



da 

+ = dz 

The first equation (|2.16p at a > has the following solution: 



a{t) = 



aa = a(0) > 0, 



(2.16) 



(2.17) 



1 + a^t 

whereas the second equation of (|2.16p corresponds to the hydrostatic pressure distribution 
in the whole space, and the sign (+) corresponds to an ascending flow (the axis z is 
directed upwards), while the sign (— ) corresponds to a descending flow (the axis z is 
directed downwards). In this case, equation (|2.12p in the region r ^ Rq has a solution 



uj{t) = w(0) exp 



ao dt' 
1 -t- aQt' 



uiO){l + aot), 



(2.18) 



which corresponds to the linear growth of vortex rotation velocity in time. 

At last, we notice that equation (|2.12p under the condition of a{t) = w(<) > takes 
the form: 

doj 



dt 



''{t)=0. 



(2.19) 



The solution of the nonlinear equation (|2.19p corresponds to the so-called "explosive" 
instability: 

u{0) 



uj{t) = a{t) 



(2.20) 



1 - w(0)t' 

when in a finite time interval t^ = 1/ijj(0) = l/a(0) the angular velocity of a fluid rotation 
uj{t) and the derivative of the axial velocity with respect to z, i.e. a{t) = dvz/ dz^ formally 
approach infinity, although they are actually limited from the above due to the effects of 
compressibility of a liquid. 

The positive sign of a corresponds to the growth of velocity of the flow along the axis 
z directed by velocity Vzq. Thus, in a descending flow of a liquid which flows out through 
a hole at the bottom under the action of gravitation and accelerates along the axis z 
by the linear law (|2.9p . we get an exponential (j2.13p . linear ()2.18p . or "explosive" (|2.20p 
laws of the acceleration of vortex rotation in time. At this, the amount of a liquid, which 
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flows out, is completely compensated by the inflow of the same amount of the liquid 
with velocity of the converging radial flow ()2.10|) from the surrounding region, which is 
considered as a large enough reservoir of the substance. 

Such a simple model explains the funnel formation in a bath at the opening of a hole 
or whirlpool formation on a river in the place of a sharp deepening of the bottom. The 
deceleration of the rate of rotation velocity growth of a liquid (water) in the stationary 
regime is caused by the friction with fixed solid surfaces as well as by the energy dissipa- 
tion on the tangential jump of the azimuthal velocity at the point r — Rq (see section 

This model may also explain the origination of sandy tornados in deserts. Due to a 
strong heating of some sites of a surface of a sandy ground by the sunlight (the darkest 
or located perpendicularly to the solar rays) , the nearby air gets warm locally and starts 
rising upwards with acceleration under the action of the Archimedean force as less heavy. 
If this acceleration in quasi-stationary conditions corresponds approximately to the linear 
in z law (|2.9p . then at a > we again obtain the exponential law of rotation velocity 
growth of a vortex (|2.13p . At this, the accelerated decrease in time of the pressure near 
the vortex axis r — Q with the increase of uj{t) and a{t) leads to the suction of the sand 
deep into the vortex and to formation of visible tornados, which are frequently observed 
in deserts. 



3. Vortices in mult i- component open systems with chemical and 
phase transformations 

Consider a multi-component system in which chemical processes or phase transitions 
are taking place. Partial balance equations of substances for each component (phase) 



look Hke (see Sedov (1997)): 



^-t-div(p,v,) = g., ^^J^QaV, (3.1) 

where pi and are density and hydrodynamical velocity of the i-th component, Mi is 
the mass of this component in a volume V, Qi is the capacity of a source {Qi > 0) or 
sink [Qi < 0) of this component due to chemical reactions or phase transformations. 

For a closed system the conservation law of the total mass of all N components holds 
true: 

N N 

M = M, = const, ^ Q, = 0, (3.2) 

1=1 i=l 

and therefore the general continuity equation can be written as 

d N N 

+ div(pv) = 0, p = ^pi, pv = ^p,v,. (3.3) 

1=1 j=i 

For open systems with a bulk sink or source of the substance and its unlimited inflow 
from the environment under conditions of dynamic and chemical balance, one or several 
components can appear in the system or drop out from the general collective motion 
of the matter (liquid or gas) as a result of chemical reactions or phase transformations. 
At the same time, diffusion and hydrodynamical flows support constant density and 
chemical potential in the open system, i.e. under stationary conditions it can be with a 
good approximation assumed that p = const. In this case the effective continuity equation 
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divv = Q/p^O, (3.4) 

where Q is a certain capacity of a source (Q > 0) or sink {Q < 0) of the substance in a 
system or in a some part of its volume V. 

Further we shall consider systems with a homogeneous in space bulk sink of matter in 
a some finite cylindrical region with the radius Rq- In cylindrical coordinates, equation 
(|3.4p looks like (at the absence of longitudinal flows, Vz = 0): 

divv . ^ + ^ = I -1^1/^ " -l/r, r < i?o, 
dr r { 0, r > Rq. 

From the last expression it follows that the radial velocity can be set as 

/ A / -Pr, r < Ro, 

where /3 = 1/2t > 0. We notice, that \Q\, t and /3 can be either constants, or variables 
in time, but must be homogeneous in space, i.e. the condition Vdiv v = 0, which is 
necessary and sufficient for the use of the Navier-Stokes equations (j2.1|) and ()2.2p for an 
incompressible liquid, is satisfied. 

Substituting expressions (|2.7p and (|3.6p for azimuthal and radial velocity profiles in 
equations (|2.ip and (|2.2p . we get the equation for determination of hydrodynamic pressure 
(compare with (|2.14p ): 



.3 V"'+^'+rf^]f ^^'^^ 



(3.7) 



and also the equation for the angular velocity of rotation of the liquid (gas) in a vortex 
(compare with (|2.12p ): 

duj _\ 2Puj, r ^ Ro, 



dt \ 0, r>Ro. ^^-^^ 

In the case of a constant sink \Q\ = const and /? = const wc get from p.Sp the 
exponential law of acceleration of a vortical rotation (at uj{0) ^ 0): 

u;(t) = w(0) e*/^ (3.9) 

If the capacity of the surrounding reservoir is finite, and the velocity of inflow (and also 
the sink) of a substance in the dynamic equilibrium state decreases in time, for instance, 
by the exponential law 

^(i)^ Iffl^^^e-*/*", i.e. r(t)=roe*/*o, (3.10) 
2 p 

then the time dependence of the angular velocity of a vortex rotation, according to p.Sp . 
looks like 

'to 



w(0=w(0)exp|-^(l-e-*/*°U (3.11) 



To 

and tends to the maximum value Wqo = "^(0) e*"/"^" at t — > oo, whereas during the time 
t <S^to the exponential growth of velocity Lij{t) « w(0) e*/"^" remains. 

From equation p.Sp it follows, that at /3 > and either sign of w, there exists another 
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non-stationary solution at /? = i|cLi(i)|, when equation (|3.8p becomes (for r ^ Ro)'- 

^--^W=0. (3.12) 

The solution to this equation corresponds to the nonlinear instability of the " explosive" 
type (compare with (|2.19p '): 

1 - \uj{0)\t 

For a finite interval of time to — l/|i^(0)|, the angular velocity w(t) as well as /3(i) formally 
reach their infinite values. 

Thus, the presence of a bulk sink together with the unlimited inflow of a substance in 
open systems with chemical reactions or phase transformations, due to nonlinear hydro- 
dynamical forces, which arise under the action of converging flows Vr < 0, results in the 
acceleration of vortices even at the absence of axial flows {vz = 0). An example of such 
a natural mechanism of vortex acceleration due to chemical reactions can be origination 
of the so-called " fiery tornados" , which arise during big fires in closed volumes at the 
formation, on the one hand, of solid combustion materials, in particular solid compounds 
of oxygen (oxides) that completely exclude oxygen from the air, and, on the other hand, 
at the presence of free infiow (draught), which creates radial and ascending flows of the 
external air enriched with oxygen. 

Moreover, during the dissolution of solid crystals in a liquid (for example, manganic- 
sour potassium in water), when the heavier solution is immersed downwards by gravity, 
and a new pure solvent with origination of converging radial flows goes to its place, a 
hydrodynamical vortex may arise, which due to friction will force the dissolving crystals 
to rotate. 



4. Origination of tornados and typhoons during formation of dense 
cloud systems 

Let us consider one more concrete example of origination of vortices in an open non- 
equilibrium heterogeneous system with phase transformations, namely, in a humid at- 
mosphere during water vapor condensation in rain clouds where tornados and typhoons 
might originate. 

Despite the colossal quantity of available observable data on tornados and typhoons in 
natural conditions and numerous attempts of modeling these powerful atmospheric vor- 
tices in laboratory plants and computer simulations (see Bengsson & Lighthill (1982)), 
the true reasons for origination and development of such phenomena are not yet un- 
derstood completely. The majority of theoretical models of tornados and typhoons can 



be reduced to the well-known "funnel effect" in hydrodynamics (see Kundu (1990) 



Sedov (1997) ). At this, it is supposed that the radial squeeze of twirled descending flows 



of dense humid air is being realized by the air masses, which are flowing in the region of 
low atmospheric pressure. 

In Ref. (Pashitskii (2002) ) an essentially new mechanism of tornado and typhoon orig- 
ination during formation of dense cloud systems was introduced. This the mechanism is 
directly connected to the process of intensive condensation of water vapor at cooling of 
the humid air below the dew-point. A "vapor - liquid" phase transition with the for- 
mation of weighed water droplets in a cloud (fog) is accompanied by a large decrease 
(about 800 times) of the specific volume which is occupied by water molecules. This 
means that the condensation of vapor should result in essential (almost twice at initial 
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100% humidity of the air) decrease of the density of a gaseous component of the two- 
phase heterogeneous system "air - water drops" in comparison with the initial density 
of a chemically homogeneous gas mixture " air - water vapor" . It is obvious that this is 
equivalent to the existence of a bulk sink (convergence) of a substance inside the cloud 
and origination of the concentration gradient of water molecules on the border of a cloud. 

As it follows from the Boltzmann equation for a mixture of gases with close by weight 
molecules (see Lifshits & Pitaevskii (1981)), the diffusion flow of one of the components 



under the action of the corresponding concentration gradient, as a result of intermolec- 
ular collisions, creates an average hydrodynamical flow of the gas mixture as a whole. 
Since the cloud is a thermodynamically non-equilibrium open system with the unlimited 
inflow of a substance from the outside, the presence of a bulk sink and the humidity gra- 
dient caused by the condensation of moisture should result in appearance of converging 
hydrodynamical flows of the humid air from the surrounding areas of the atmosphere. 
Under conditions of a dynamical balance, such flows retain a constant in time and almost 
homogeneous in space average density p — const, which is equivalent to the condition of 
incompressibility of the medium. 

As has been shown in the previous section, within the framework of the hydrodynamical 
model of an incompressible viscous liquid with a bulk sink (convergence) of the substance 
and converging radial flows, under the joint action of convective and Coriolis forces an 
instability of a vortical "rigid-body" rotation appears, which can evolve in time either 
by the exponential law (in the case of constant sink and inflow of substance), or by a 
scenario of the " explosive" type when the infinite rotation velocity is reached in a finite 
time interval (under the condition of simultaneous unlimited increase of velocities of the 
sink and the inflow). 

We will show, that such an instability can be one of the reasons of origination and 
acceleration of powerful tornadoes and typhoons during the formation of dense cloud 



systems (see also Pashitskii (2002) I 



4.1. Balance of substance during the condensation of rain clouds 

The rain cloud represents an open non-equilibrium two-phase system (a mixture of the 
air and a fog) that self-organizes due to condensation of water vapor and keeps small 
water drops in a certain volume by converging flows of the air and the Stokes friction 
forces. 

Consider a cloud of a cylindrical form of the radius Rq and height h. As inside a 
forming cloud the process of vapor condensation from the humid air, cooled below the 
dew-point Tjcw (at the given pressure distribution P) takes place with the certain ve- 
locity, the effective balance equation of the substance inside a cloud, with account of 
horizontal (radial) and vertical (ascending or descending) air flows from the surrounding 
atmosphere, and a bulk sink of the substance with capacity Q < 0, can be written as 

2TTRohvr{Ro) + TTRl[v^{h) - i;^(0)] = ^^-^TiRlh. (4.1) 

P 

Substituting expressions (|2.9|) and (|2.10p in equation (|4.1I) . we get the relation 

2/3 = a+M, (4.2) 
P 

which corresponds to the effective continuity equation (compare with (|2.4p and p.Sp ) 
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Thus, at a > there is a strengthening of the nonhnear instabiUty, which results in 
acceleration of a vortex rotation in the region r ^ i?o by the following law: 



uj{t) = w(0) exp 



dt' 



1 



Tit') 



(4.4) 



At the same time, in the external region r > Rq, in agreement with (|3.8p . the angular 
velocity remains constant and is equal to the initial value a;(0). This means that on the 
border of the vortex core at r = i?o a tangential discontinuity of azimuthal velocity arises 
and grows in time. This jump at a = const and |Q| = const equals 



Av^{Ro,t) = Voit)^u;{0)Ro[e^''' 
and in the case of an " explosive" instability we have 

1 



1], 



Vo{t) ^ ioiO)Ro 



1 - Lu{0)t 



1 



(4.5) 



(4.6) 



4.2. Instability of the tangential jump of velocity on the border of the tornado core 



As is well-known (see Landau & Lifshits (1987) I, a tangential jump of velocity in hydro- 



dynamics results in origination of instability of surface perturbations of an incompressible 
fluid. The dispersion equation of small perturbations in linear approximation with ac- 



count of a finite viscosity looks like (see Pashitskii (2002) I : 

{uj — kVo){uj — kVo + ivk^) — —uj{uj + ivk"^), 



(4.7) 



where k = 27r/A is the wave number of surface waves with the length A, which are 
distributed along a flat surface. The equation (|4.7[) with a good approximation can be 
also applied to the case of the cylindrical surface of a radius i?o if we consider A;i?o ^ 1- 
The frequency ujQ{k) = Re uj{k) and the increment 7(fc) = Im w(fc) of such short-wave 
surface perturbations on the tangential jump of velocities with transversal displacements 
of a surface along the z axis and the propagation along the x axis (at fc > 0) 



C(a;, z, t) = (o e e'^* exp(iA:a; - iuJot), 
according to dispersion equation (|4.7p . are equal to 



^o(fc) = — , lik) ^ - 



kW^ + v^k^- vk' 



(4.8) 



(4.9) 



If fc < V^/v (but fci?o > 1), then the increment is 7(fc) = kVo/'i. At fc > Vq/v, 
according to (|4.9p . the increment tends to its maximum value (at k —>■ oo): 



7o 



(4.10) 



It is worth noting, that relations (|4.7p - (|4.10p were obtained under the assumption of a 
constant value of the velocity jump Vb = const. But if the condition ^ 2/3 is satisfied, 
i.e. if the jump of the velocity (j4.5p grows in time very slowly (adiabatically) in comparison 
with development times of the unstable short-wave excitations with increment (|4.10p . it 
is possible to present approximately the time evolution of the amplitudes of surface 
perturbations at z = in the form 



IC(i)I^Coe 



Co exp 



Lu^{0)Rlt 
2^' 



2l3t 



-'1 



(4.11) 
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for the exponential instability of a vortex (for f3 — const), or 



|C(i)| =Co exp. 



u\0)Rlt 
2v 



1 

1 



1 - w(0)i 



2 



(4.12) 



in the case of the "explosive" instability given by (12.191) and (I3.13p . 
During the time interval t <C 1/2/3, when Vo{t) ~ 2(3uj{0)Rot, we have 

mi ~ Co exp , (4.13) 



instead of the "exponent to the exponent" law (|4.1ip . 

A more detailed analysis of time evolution of small surface perturbations at a variable 
velocity Vo{t) can be obtained in the case where v = with the help of the following 
linear differential equation: 

S - fc^o W§ + ^fc'K,' WCW = 0. (4.14) 

If Vo(t) ~ i, then the dominate time dependence of the solution (|4.14[) looks like 

C(t) exp {at^ - fot} , (4.15) 

i.e. it is slower, than the dependence (|4.13p . But if one takes into the account the fact, 
that at ^ the inequality k <^ Vo{t)/i' always takes place, i.e. the increment is 
7 = kVoit)/2, then instead of ([4A3ll we get (at 2/3t < 1): 

\C{t)\^Coexp{f3u{0)kRot^}, (4.16) 

i.e. the use of the "adiabatic expressions" (|4.1ip is possible at early stages of a vortex 
instability evolution. Similarly, in the case of "explosive" instability, with the account of 
(|4.6p and (I4.12|) under the condition oi oj{Q)t ^ 1, we obtain 

\at)\-Coe^p{kLo\0)Rot^}. (4.17) 

4.3. Development of the turbulence on the border of the tornado core 

As is known, the maximum value of the growth increment of unstable perturbations in 
a some non- homogeneous layer of thickness I corresponds to the wave number k^ ~ 1//. 
In the case of surface waves on a tangential jump of velocity, the role of the thickness of 
a transition non-homogeneous layer is played by the doubled amplitude of these waves 
(surface displacement), i.e. km — l/2\((t)\. Substituting this value in expression (|4.16p . 
we get the equation 

InKM^MW^ (4.18) 

Co mt)\ ^ ' 

from which it follows that \C{t)\ changes approximately according to the law \C{t)\ ~ 
t'^/lnt. At this, the effective coefficient of the anomalous turbulent viscosity inside a 
surface layer of thickness 2|C(i)| can be estimated as 

,,*{t) C\t)jit) ^f3u;{0)Rom\t, (4.19) 

i.e. the value iy*{t) grows in time, according to (|4.18p . almost by the cubic law (^ t^) 
and can reach very large values {v* :s> v). 

The saturation of turbulent perturbations occurs due to the nonlinear dissipation, 
when the increase of the kinetic energy of surface waves, as a result of instability of the 
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tangential jump of velocity ^jpv^, is compensated by the turbulent dissipation energy 



2 

per unit time and unit volume pv^ /I (see Landau & Lifshits (1987) I, where v ~ d\C,{t)\/dt 



is the velocity of turbulent pulsations, I ~ |C(^)I is a characteristic scale of turbulence, 
and 7 ~ 7rVo/|C| is the maximum increment of the instability. Therefore, it follows that 

vit) c ^Vo{t). (4.20) 

On the other hand, the saturation of acceleration of the angular velocity of a "rigid- 
body" rotation of the vortex core (tornado) occurs when the azimuthal velocity v,p(t) = 
RoUj{t) and the value of a velocity jump Vb(t) reach velocities close to the speed of sound 
in the air Cg — 330 m/s, i.e. when the compressibility effects of the air and the effects 
of a finite bulk viscosity appear. The characteristic time of a rotation saturation of the 
vortex, according to (|4.5|) , equals 



w(0)i?o 



(4.21) 



Thus, Wniax — Cs and 7max - TTCs/Cmax, where the value of (max, according to l|4.18p . is 
determined by the equation: 

MO^Ro^ (4.22) 

Co 2Cniax 

Let us notice that in the case of the "explosive" regime of acceleration of a vortex, it 
is necessary to replace /3 by lo{Q) in expressions (|4.18[) . (|4.19[) and (|4.22[) . and the value 
^max under the condition ^ aj(0)i?o equals t^ax — 

If we assume that the initial azimuthal and radial velocities of air on the borders of 
a cloud, which has radius i?o — 1 kui, are equal by order to = w(0)i?o = PRo — 10 
m/s, i.e. /3 ~ 10^^ s^^, then the time of acceleration of the tornado core to the speed 
of sound Cs equals tmax — (1-^2)- 10^ s for the "explosive" or exponential acceleration 
laws of a tornado, respectively. Setting the minimum initial value of displacement Co of a 
surface r = i?o to the free path of molecules in the air at the normal pressure — 5 • 10~^ 
cm, we get from (j4.22l) the following estimation for the maximum scale of the turbulent 
pulsations: Cmax — (20 80) m. 

In this case the maximum value of the turbulent viscosity in the surface layer of a 
thickness 2Cmax on the border of a vortex core r = Rq equals 

^max - ^^^(0)i?oCmaximax/3 (4.23) 

and can by many orders exceed the usual coefficient of the kinematic viscosity of the 
air. In particular, for the above mentioned estimations imax and Cmax we get the value 
j^j^jjx — 5 • 10^ cm^/s, whereas the usual kinematic viscosity of the air equals ly ~ 0.15 
cm^/s, i.e. their ratio is v^^^^/v — 3 • 10''. 

Thus, on the surface of the central part of a powerful tornado of a radius i?0 7 which 
rotates almost with velocity of sound V^max = ^{^)Ro exp[7niaximax] — Cg, a very viscous 
layer of a thickness 2Cmax ^ and with the viscosity v^^^ » v appears because of 
the turbulence. Due to this fact, a partial entrainment of the nearby air in the region 
r > Rq occurs, and therefore the profile of azimuthal velocity distribution of the air in 
this region changes. We will seek it in the form 

v^{r,t)=n{r,t)Rl/r {r > R), (4.24) 



where n{r, t) is an unknown function of r and t, which, according to (|2.2p . satisfies the 
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equation: 



an (3R^dn _ , fd^n ion 



dt r dr '"^'^ V 9r2 r dr 



^maxUlT--7r- (4.25) 



with a boundary condition f7(i?o,t) = uo{t) — Lo{i)) e^''*. From here we get the following 
expression for the rotation velocity of the air in the external region r > Rq : 

.,.(.<,=.(0,«.e-(^)""«I, (4,20) 

where Ka{x) is the McDonald function (the modified Bessel function), a = {Ro/'2,l*Y — 1, 
and I* = y/ v^^y./2f3. Using the obtained above value of i^^ax; get the following 
estimation: I* ~ 150 m. 

On large distances r ^ I* the velocity (|4.26p changes according to the law 

vAr) - [^y" ^ e-^f^\ (4.27) 

that explains rather a weak motion of the air around a tornado on the distances of r > 150 
m. On the other hand, those fact that the turbulent pulsations and rotation velocity of 
a vortex core can reach the velocity of sound ("Dmax — V'^max — Cs) explains the reason of 
the intensive generation of low-frequency sound waves in powerful tornadoes (a "roar" 
of a tornado) . 

4.4. Formation of a tornado funnel with account of gravity and vertical flows of the air 

As was mentioned above, for the flows of the air, which flow into a rain cloud of a cylin- 
drical form during the time of its condensation, the following hydrodynamical velocities 
are natural: 

-/3r, r ^ Rq, _ f wr, r Rq, _ ( w^o + az, r ^ Rq, 

-pRl/r, r>RQ, '"'^ ' \ uR^r, r > Rq, ' \ 0, r > Rq, 

(4.28) 

where the parameter (3 is determined by relation (j4.2p . In the case of exponential insta- 
bility we have Lj{t) = uj{0)e'^'^* (for a = const and \Q\ = const). At this, bulk viscous 
forces in equations (|2.ip - (|2.3p are equal identically to zero, so that we get 

gp ( pr[L0^it)-P^], r^RQ, 

^K(0)+/3^], r>Ro, ^'-'^^ 

9P ^ f -pg-pa^z, g^g + av^o, ri^ Rq, 

dz \ -pg, r> Rq. ^ ■ ' 

Integration of equations (|4.29p and (|4.30p determines the difference of the pressure 
between two arbitrary points: 

-p~g{z, ~zi)~^ o?[zl - ^?) + f {^\t) ~ P'] {rl - rj), Rq, 

= ^ -pg(^, - zi) - ^ [u^'iO) + /32] [^~^)^ ^1,2 > Ro- 

(4.31) 

From (|4.3ip it follows, that the form of a surface of constant pressure (an isobar), which 
corresponds to the point of water drops evaporation Pgvp in the internal region r ^ Rq 
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is determined by the equation 

z'{r,t) + %z(r,t) + 4 [-^(t) +-^(0)] - 4 [-'W " ^ "/-^^ = 0, 

(4.32) 

where the coordinate z is counted from the initial flat surface Pq = Poo = Pcvp, smd in 
the external region r > Rq is determined by the relation 

. - m±fl^ + (4,33) 

where the function ZQ{t) is to be obtained from the condition of the isobar continuity 
at the point r — Ro. According to equation (I4.32p . we find the value of the coordinate 
z{r, t) at the point r = Rq : 

- const, (4.34) 




^ + p^] - ^""P^ = const, (4.35) 

2.9 pg 

i.e. in the external region r > Rq^ the form of the isobar does not depend on time. 

The coordinate of a point of the isobar on the vortex axis r ~ 0, according to (|4.32[) - 
(|4.35p . is determined by the expression 




.(0, t)^^JL + ,l^- ^°[-^(^);-^(")] + 2(^_^. (4.36) 

From the last formula it follows, that for 5 > 0, with the growth of the angular velocity 
Loify as a result of the exponential instability, there is an increase of the absolute value 
of the negative coordinate 2;(0,t), which corresponds to the deepening of the minimum 
of the function z{r^ t) in the region z < 0. 

In Fig. [TJ cylindrically-symmetric surfaces z{r, t) — z{r, t) — z(oo, 0) and their sections 
by mutually perpendicular planes depending on r for the isobar P = Pcvp, which corre- 
sponds to the bottom edge of a condensation area of the moisture at r — > cx), are shown 
for 5 consecutive moments of time with the identical intervals. They display the time 
dynamics of a funnel evolution, which is filled with a fog (a heterogeneous mixture "air 
— water drops"), and develops during the origination of a tornado on the bottom edge 
of a cloud. The deepening of the funnel occurs till the contact of its minimum that cor- 
responds to the value Zmin — —g/o? with the surface of the Earth, or till that moment 
of time, when the maximum velocity of a vortex rotation reaches the speed of sound. 

Thus, the account for the gravity and vertical flows allows to describe one of the main 
observable phenomenon — a funnel formation on the bottom edge of a cloud during the 
origination and development of tornados. 

4.5. Discussion of viability of the obtained results for the description of tornados and 

typhoons 

With the purpose of establishing of the viability of the obtained new non-stationary 
solutions of the Navier-Stokes and continuity equations for an incompressible viscous 
medium with a bulk sink and free inflow of the substance for the description of atmo- 
spheric vortices (tornados and typhoons), we will carry out numerical estimations of the 
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Figure 1 . Dynamics of the isobar shape on the bottom edge of the condensation area during 
formation of a cloud and development of the hydrodynamical instability of a vortex (tornado) . 
Different (descending) shapes correspond to different equidistant times of vortex evolution. 



characteristic times of development of hydrodynamical instabilities of vortical motion 
under conditions of intensive condensation of moisture inside a cloud. 

Let us consider a round cylindrical cloud of a radius Rq ~ 1 10 km slowly rotating in 
the twirled air flow with equal by order initial values of azimuthal and radial velocities on 
the border of a cloud Vro ~ Wi^o ~ (1 ^ 10) m/s. This corresponds to the absolute values 
of initial angular velocity a;(0) and velocities of a converging radial flow /3 in the range 
(10~^ 10~^) s~^. For humidity of atmospheric air about 100%, its density almost twice 
exceeds the density of the dry air po = 1.3 x 10"'^ g/cm^, so that the mentioned values 
of the parameter /? are equivalent to capacities of a bulk sink Q due to condensation 
of a moisture of the order of ~ 5 x (10"'' -t- 10~^) g/s-cm^. Under these conditions, the 
characteristic time of acceleration of a vortex rotation in the case of the exponential 
instability equals t = 1/2/3 « (1 100) min. 

The maximum velocity of the order of Cg is reached in a time interval of t » (3 500) 
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min. Typical observable times of origination and development of a tornado indeed lie in 
this time interval (from several minutes to several hours). 

For large-scale atmospheric vortices such as cyclones, hurricanes and typhoons, which 
originate in cloud masses with sizes of about 100 km and more, characteristic times of 
the instability development increase in two - three orders of magnitude and can last 
for several days, which also agrees with the observable times of hurricane and typhoon 
existence. At the same time, the decreasing pressure on the axis of a vortex caused by the 
cyclostrophic rotation regime explains the characteristic "sucking" effect of a tornado. 

Moreover, as marked above, in a two-phase system "air — water drops" the decrease 
of the pressure below the boiling point of water at the given temperature should result in 
the termination of the process of moisture condensation and lead to evaporation of the 
droplets. This explains the formation of a clear tornado "core" or a typhoon "eye" in the 
central part of a cloud system in the paraxial area of a vortex. Inside this area the velocity 
of a vortical rotation of the air slows down in time, since at the condition of P < Pcvp 
there is a bulk source of the gas phase {Q > 0) instead of a bulk sink (Q < 0), and the 
parameter /3 changes its sign (/3 < 0), that corresponds to the exponential deceleration 
of the air, according to (|3.8p . Such a condition of a dead calm is observed in the center 
of the typhoon " eye" . 

It is also possible to estimate the moment of time tn when a tornado funnel contacts 
with the surface of the Earth. In the case of the exponential instability, at z = —H 
(where H is the height of the bottom edge of a cloud above the ground), according to 
equations (I4.36P and ()3.8|) . under the conditions oj'^it) :§> and g/a'^ ^ H we find 



iff « — In 

4/3 



~9H 



(4.37) 



Ai H — 1 km for the above mentioned values of v^q we get the estimation oltn ~ (2-^400) 
minutes. 

In conclusion of this section it is necessary to emphasize, that the growth of the kinetic 
energy and angular momentum of a vortex during the development of the exponential 
or "explosive" instability does not contradict to the conservation laws of energy and 
angular momentum, because in an open system with a bulk sink and continuous inflow 
of the substance, the transfer of the necessary amount of the kinetic energy and angular 
momentum to the region r ^ i?max from the environment with a nonzero initial rotation 
velocity v^q = aj(0)i?g/r takes place. The value of the maximum radius i?rnax in natural 
conditions at tornado and typhoon origination on a free space of the land or sea is 
actually limited only by a finite curvature of the Earth's surface. This means that the 
maximum magnitudes of energy and momentum of the atmospheric vortex cannot exceed 
the following values (per unit length): 

EMn = np / rdr [«2 (r) -1- v^^{r)] = npR^ [lj^{0) + P^] In (i?,„ax/i?o), (4.38) 
•I Ro 

M,_ = 27rp / r^v^dr = i: pRlRl^^^uj'^ {Q) . (4.39) 

For the above mentioned values of the parameters typical for tornados and for -Rmax = 
10^ km we get E^^y, « (10*^10^2) At this, the maximum intensity liberated during 

the acceleration of a vortex \ ^ 4/3£'max reaches roughly 1 GW per meter of 

V / max 

the vortex length, which explains the enormous destructive force of a tornado. 

Thus, the carried out numerical estimations show that the introduced in the paper a 
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new mechanism of hydrodynamical instability of a vortex in an incompressible viscous 
medium under the action of convective and Coriolis forces created by converging radial 
flows in open thermodynamically non-equilibrium systems with a bulk sink and unlimited 
inflow of a substance from the environment may be the real cause of origination and 
development of powerful atmospheric vortices — tornados and typhoons — during the 
intensive condensation of water vapor from the humid air cooled below the dew-point 
during formation of dense clouds. 



5. Conclusions 

In the present paper a new class of exact solutions of hydrodynamic equations for 
an incompressible liquid (gas) at the presence of a bulk sink and ascending flows of a 
substance was considered. For existence of such solutions it is essential that one or several 
components of a multicomponent fluid are excluded from the collective hydrodynamic 
motion of the liquid at the expense of chemical or phase transitions, but due to dynamical 
and chemical equilibrium in the open system with the surrounding medium, a constant 
in time and almost homogeneous in space chemical composition of the matter as well as 
constant density p — const are maintained. It is shown that those profiles, which nullify 
the terms in the Navier-Stokes equations that describe viscous effects, exist and represent 
vortex structures with "rigid-body" rotation of the core and converging radial flows. In 
the case of constant bulk sink and inflow of the matter from the outside, the azimuthal 
velocity of a "rigid-body" rotation increases exponentially in time. At simultaneous 
infinite increase of the sink and inflow rates, Vip increases by a scenario of the "explosive" 
instability, where during a finite time interval an infinite value of rotation velocity is 
reached. 

On the basis of the developed theory of unstable hydrodynamical vortices, the offered 
in Pashitskii (2002) mechanism of origination and development of powerful atmospheric 
vortices — tornados and typhoons — during the intensive condensation of water vapor 
from the cooled below a dew point humid air during the formation of dense rain clouds 
is considered. Within the framework of this mechanism it is possible to explain the basic 
characteristics of tornados and typhoons. The characteristic times of development of 
the instability of a vortical motion agree in order of magnitude with the corresponding 
times of origination and existence of tornados (from several minutes to several hours) 
and typhoons (several days). The acceleration of a vortex rotation is limited by velocities 
comparable with the speed of sound, when effects of the compressibility of air, dissipative 
forces of a bulk viscosity and large-scale turbulence start playing their roles. 

With the account of gravity, the proposed model describes the main phenomenon of 
a tornado — the formation of a lengthening funnel on the bottom edge of a cloud as a 
result of the change of the shape of the surface of a constant pressure (isobar) that bounds 
from the below the area of intensive condensation of moisture. It is worth noting, that 
in an open system there are no problems with the conservation laws of the energy and 
angular momentum, because the radial flows transfer not only the necessary amount of a 
substance, but also the necessary amount of the kinetic energy and angular momentum 
from the slowly rotating environment. 

The velocity jump on the border of a tornado core results in the development of a 
strong turbulence, which can be described with the help of the anomalous coefficient 
of turbulent viscosity, which in many orders exceeds the usual viscosity of the air. Our 
estimations of the turbulent viscosity coefficient are agreed with the values known from 



the literature (see e.g. Kochin, Kibel & Rose (1964)). 



18 



E. A. Pashitskii, V. N. Malnev, and R. A. Naryshkin 



REFERENCES 

Nalivkin. D. V. 1969 Hurricanes, Storms and Tornadoes. Leningrad, Nauka (in Russian). 

Bengssok, L. & LiGHTHiLL, J. (eds.) 1982 Intense Atmospheric Vortices. Berlin. 

KuNDU, p. K. 1990 Fluid Mechanics. Academic Press Inc. 

Landau, L. D. & Lifshits, E. M. 1987 Fluid Mechanics. Pergamon Press. 

Sedov, L.I. 1997 Mechanics Of Continuous Media. In 2 vol. World Scientific. 

GOLDSHTIK, M. A., SiiTERN, V. N. & Yavorsky, N.I. 1989 Viscous Flows with Paradoxical 

Properties. Novosibirsk, Nauka (in Russian). 
Pashitskii, E. A. 2002 Appl. Hydrodyn., 4, 50 (in Russian). 
Lifshits, E.M. & Pitaevskii, L.P. 1981 Physical Kinetics. Pergamon Press. 
KoCHiN, N.E., KiBEL, LA. & Rose, N.V. 1964 Theoretical Hydromechanics. Wiley- 

Interscience. 



